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ABSTRACT 

We have observed four X-ray underluminous groups of galaxies using the 
Giant Meterwave RadioTelescope. The groups NGC 524, 720, 3607, and 4697 
are underluminous in relation to the extrapolation of the L x - T relation from rich 
clusters and do not show any evidence of current AGN activities that can account 
for such a departure. The GMRT observations carried out at low frequencies (235 
and 610 MHz) were aimed at detecting low surface brightness, steep-spectrum 
sources indicative of past AGN activities in these groups. No such radio emissions 
were detected in any of these four groups. The corresponding upper limits on 
the total energy in relativistic particles is ~ 3 x 10 57 erg. This value is more 
than a factor of 100 less than that required to account for the decreased X-ray 
luminosities (or, enhanced entropies) of these four groups in the AGN-heating 
scenario. Alternatively, the AGN activity must have ceased ~ 4 Gyr ago, allowing 
the relativistic particles to diffuse out to such a large extent (~ 250 kpc) that 
their radio emission could have been undetected by the current observations. If 
the latter scenario is correct, the ICM was pre-heated before the assembly of 
galaxy clusters. 

Subject headings: galaxies : clusters — clusters : individual (NGC 524, 720, 
3607, 4697) — galaxies : active — galaxies : intergalactic medium 



1. INTRODUCTION 

Recent studies of the hot intracluster medium (ICM) in clusters of galaxies have shown 
that this gas has been heated by some non-gravitational sources other than the heating ex- 
pected from falling into the gravitational potential of dark matter. In particular, the entropy 
of the gas, defined as S = T/nl^ 3 , where T and n e denote the gas temperature and the elec- 
tron number density respectively, appears to be in excess of expectations from gravitational 
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interactions of gas with dark matter. This excess entropy makes the gas less luminous than 
expected and affects a deviation from the scaling relation between X-ray luminosity and 
gas temperature, especially in poor clusters (Ponman et al. 2003, and references therein). 
There have been a number of proposals to explain this entropy floor including energy input 
from supernovae (e.g., Wu et al 2000), warm-hot intergalactic medium (Valageas et al 2003), 
radiative cooling (e.g., Voit & Bryan 2001), accretion shocks (e.g., Tozzi & Norman 2001), 
and AGN heating (e.g., Nath & Roychowdhury 2002). 

One of the theoretical models to explain the X-ray observations invokes energy input 
from active galactic nuclei (AGN). It has been proposed that AGN heating not only can slow 
down or even balance the precipitous cooling of the gas in the central region, but it can also 
raise the gas entropy outside the core (r > 100 kpc). These models were partly motivated by 
the presence of X-ray-deficient bubbles in the inner regions of many clusters, for example, the 
Hydra A cluster (McNamara et al. 2000), or in the periphery of some clusters, for example, 
in MKW3s (Mazzotta et al. 2002). 

Such bubbles, presumably fed by AGN activity, are expected to rise in the cluster 
atmosphere because of buoyancy and as they expand in size, they are expected to do work 
against the ambient gas, providing a source of heating for this gas. Begelman (2001) and 
Ruszkowski & Begelman (2002) suggested this 'effervescent heating' could explain the related 
problem of quenching cooling flows in galaxy clusters. Roychowdhury et al (2004, 2005) 
studied the effect of 'effervescent heating' outside the core and found that it could explain 
the X-ray observations. According to their findings, a cluster with a total mass ~ 10 14 M , 
with gas at ~ 1.5 keV, requires an AGN, or a number of AGNs, that could supply energy of 
order 5 x 10 61 erg over a time scale of ~ 1 Gyr (Roychowdhury et al 2005, Fig. 4). 

Recently, Croston et al (2005) studied the X-ray luminosity and temperature relation 
in a sample of poor clusters after having divided the sample into two groups : clusters with 
and without radio-loud objects in them. They found tentative evidence that clusters with 
radio-loud objects in them deviated more from the scaling relation than clusters without 
radio- loud objects (Fig. 1). This lent credence to the idea that AGNs are responsible for 
the excess entropy in poor clusters. Croston et al (2005) estimated the energy required to 
lower the entropy of a group to sufficiently change its position in the X-ray luminosity and 
temperature plot, to be > 10 58 ~ 59 erg. 

However, if one extrapolated the cluster L^-T relation, it would fall above all data points, 
indicating that all poor clusters in their sample deviated from the scaling relation, although 
clusters with radio-loud objects in them deviated the most (Figure 1). Such a deviation is 
clearly indicative of factors other than only the gravitational interactions responsible for the 
entropy of the ICM gas (Roychowdhury & Nath 2003). 
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Theoretical modeling of the AGN heating by Roychowdhury et al (2004, 2005) suggests 
that the gas entropy can remain high for a considerable period of time (t ~ 10 10 y) even 
after the AGN activity has ceased. It is therefore reasonable to assume that the reason for 
the increased entropy in some galaxy groups can be the past activity of an AGN in them, 
and not necessarily the current phase of activity. One way to probe the past activity of an 
AGN is to study the low energy relativistic particles via synchrotron radio emission at low 
frequencies (< 1 GHz). The synchrotron life time of relativistic electrons with ~ 10 GeV, 
which can radiate at ~ 600 MHz for B ~ 1 //G, is ~ 1A(E/10 GeV) _1 (5/l fiG)~ 2 Gyr. 

We selected four groups of galaxies, viz. NGC 524, NGC 720, NGC 3607 and NGC 
4697, from the sample of radio-quiet groups in Croston et al (2005), which were among the 
most deviant groups from the L^-T scaling relation (Fig. 1). In the context of the AGN 
heating scenario, one way to understand the deviant radio-quiet groups is to propose the 
existence of low surface brightness and steep-spectrum radio emission in them which have 
not been detected by the existing surveys. Amongst the 4 groups mentioned, 3 have not 
been observed at low (~ 300 MHz) frequencies. It therefore appeared that these four groups 
(NGC 524, 720, 3607, 4697) were ideal targets for a deep radio continuum imaging at low 
frequencies with the Giant Meterwave Radio Telescope (GMRT). 



2. OBSERVATIONS 

The four groups NGC 524, NGC 720, NGC 3607 and NGC 4697 were observed at 235 
and 610 MHz using the Giant Meterwave Radio Telescope. Some details of the observations 
are given in Table 1. A 6 hour synthesis was carried out in each case, ensuring good visibility 
coverage so as to enable detection of anticipated weak and extended radio emission. The 
visibility coverages of these observations were adequate to detect sources up to an extent of 
~ 40'. The data analysis was carried out using the Astronomical Image Processing System 
(AIPS). On an average, ~ 10 % of the data was lost to radio frequency interference (RFI) at 
235 MHz. The data lost to RFI was much less than this at 610 MHz. The edited data was 
imaged after gain and bandpass calibration. The usual polyhedron imaging was carried out 
(using standard tasks in AIPS) to account for the well-known effects of large fields of view 
at low frequencies and the different heights of the antennas (the so-called 'w-term'). Several 
rounds of self-calibration and imaging were carried out to converge on the best images. The 
RMS achieved in each case is given in Table 1. In all cases, the RMS values in the images 
were dynamic range limited to ~ 1 part in 1000 of the peak flux density in the respective 
images. As an illustration of the quality of the images obtained, images at 235 and 610 MHz 
of a small (but, the same) region of the sky near the group NGC 4697 are shown in Figs. 2 



-4- 



& 3. 

No diffuse radio emission was detected in any of the groups to a 5er surface brightness 
limit of ~ 10 mJy/beam at 235 MHz. The 5a surface brightness detection limit at 610 MHz 
was ~ 1 mJy/beam. These images made at a resolution of ~ 20" were further smoothed 
to a resolution of ~ 1' in order to detect extended features. The smoothed images at 235 
MHz had an RMS of ~ 12 mJy arcmin" 2 . No diffuse radio emission was detected in these 
smoothed images. 



3. DISCUSSION 

We discuss the implications of our observations with regard to the energetics of rela- 
tivistic particles in the observed groups and the role of heating of diffuse gas by relativistic 
particles from AGNs. Assuming synchrotron emission at a given frequency v from relativistic 
electrons emanating from a region of length scale /, the minimum energy in relativistic parti- 
cles (electrons and protons) can be estimated from the standard minimum energy arguments 
(e.g., Longair 1994, p 292) as, 

E re ^6xl0 41 erg^ 7 LrC^ 4/7 , (1) 

where us = u/10 s MHz, L v is the power radiated by the source at v (in W Hz -1 sr _1 ), 4 pc 
is length scale in kpc and rj is the ratio of energy in protons to electrons. Converting L v 
to observed surface brightness, l v = L u (V/l) (mJy arcmin" 2 ), where V ~ I 3 is the source 
volume, this gives, 

E rel « 3.7 x 10- erg ( V ( h _)Vr (-L\ 17p (JL) 4/ \ (2) 

S V235MHz/ VmJY/arcmin 2 VIOOkpc/ V100/ K) 

It is important to recall that the standard derivation integrates the radio spectrum from a 
minimum to a maximum observed frequency, typically from 10 MHz to 10 GHz. It is also 
assumed that particles and magnetic fields fill the source volume homogeneously and the 
field is assumed to be completely tangled without any preferred direction. 

If instead of using a fixed frequency interval one had used a fixed cut-off in the energy 
of the relativistic particles, the results would be somewhat different. Beck & Krause (2004) 
re-derived the minimum energy/ equipartition energy argument by assuming a relativistic 
proton spectrum which has a spectral index 7 ~ 2.5 down to energy E p = m p c 2 , and with 
a spectral index of zero below this (but above a low cutoff at, say, Ei). This assumption 
was motivated by the observed spectrum of protons in our Galaxy. They assumed a fixed 
ratio between the number of relativistic protons and electrons, of order K ~ 100, again 
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motivated by observations in our Galaxy at a few GeV. According to them, one has for for 
the total energy density in relativistic particles, at equipartition and for a spectral index of 
7 = 2.5, 

This yields a total energy in relativistic particles (E re i = e re // 3 ), 

E„, » 1.4 x 10- erg (^-f ( h^f . (4) 

rel 6 \235MHz/ VmJY/arcmin 2 VIOOkpc/ VlOO/ V ; 

We find that these two estimates (from equations 2 and 4) differ only by factors of order 
unity. We choose to use the first equation to derive a conservative upper limit on the energy 
density of relativistic particles based on the null results of our observations. 

To derive a limit, we need to find a suitable estimate of /, the extent of the relativistic 
particles. In the absence of any other indicators, we assume it to be comparable to the 
extent of the X-ray observations, namely, the ROSAT extraction radius (r cut , as described 
in Osmond & Ponman (2004)). The values of r cut of the respective groups are given in 
Table 1. We find that for our sample, the value of r cut ~ 60 kpc. The corresponding upper 
limit on the total energy of relativistic particles for the four groups observed with GMRT 
is therefore 3 x 10 57 erg (eq. 2). Since there can be uncertainties of order unity stemming 
from a possible non-spherical geometry of the source, and other considerations, we can put 
a firm upper limit on the total energy in relativistic particles to be ~ 10 58 erg. Note that 
an T cut of 60 kpc at the mean distance to these four groups corresponds to ~ 10' and well 
within the largest angular structures (~ 40') up to which the current GMRT observations 
were sensitive. 

The groups in our sample have gas temperature in the range of 0.3-0.5 keV, implying 
^500 ~ 3.1-7.7 x 10 12 M , using the scaling relation from Finoguenov et al (2001). This 
implies virial mass M V i r ~ 10 13 M , for a median value of M500 ~ 5 x 10 12 M , assuming 
a universal profile of dark matter (Navarro, Frenk & White 1997). It has been argued that 
to reconcile with the X-ray observations of excess entropy in poor clusters, it is required to 
inject energy of order ~ 0.5-1 keV per particle (Cavaliere et al 1999; Borgani et al 2002). The 
number of particles in a M vir ~ 10 13 M cluster is ~ (il b / 'Q m ) M vir ~ 0.1AM vir ~ 1.7 x 10 69 , 
and the total energy required to elevate its entropy to the observed level is then estimated 
to be ~ 1.5-3 x 10 60 erg. 

Among the various possibilities of the basic mode of energy deposition, AGN heating 
through 'effervescent heating' has been recently studied in some detail — with bubbles of 
relativistic particles, that are created by jets of AGN when they are active, and which detach 
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from the host galaxy after the cessation of jet activity, and rise in the cluster atmosphere 
because of buoyancy, expand and do pdV work against the ambient medium (Begelman 2001; 
Roychowdhury et al 2004, 2005). Roychowdhury et al (2005) calculated the required energy 
from AGNs, which are assumed to be active for a period of 5 x 10 8 -5 x 10 9 yr, after which 
the ambient gas is allowed to evolve for a Hubble time at which point its entropy (at two 
radii, 0.1r 20 o and at r 50 o) is required to match the observed values. If the lowest envelope 
of Figure 4 in Roychowdhury et al (2005), representing the minimum energy required from 
AGNs to enhance the entropy of clusters, is extrapolated down to this mass range, then the 
energy required is ~ 10 60 erg. There is, therefore, a mismatch of a factor of ~ 100 in the 
total energy estimate and the limit from observation. 

It could be argued that the observational limit is only a lower limit, since it is based 
on equipartition, and is, therefore, close to the minimum amount of energy that relativistic 
particles can harbor. In this regard, we note that the energy estimate scales as B 2 for 
B S> B eq , where B eq represents the equipartition magnetic field strength, and as £>~ 3 / 2 
for B <C B eq . To reconcile with the above energy estimate, the magnetic field is required 
to be either B/B eq ~ 10~ 4 / 3 or B/B eq ~ 10 to reconcile with the energy mismatch. The 
equipartition value, for an energy density of 4 ~ 10~ 14 erg cm" 3 (eq. 3), is B eq ~ 1 /zG. 
Although there are uncertainties in the estimated strengths of magnetic fields in clusters, it 
is difficult to presume that the strength would differ by an order of magnitude from this value, 
as is required by the present observations (see also De Young 2006). Moreover, although the 
factor r] could be somewhat larger than the fiducial value of ~ 100 (Fabian et al 2002; De 
Young 2006), with protons dominating the energy budget, it is unlikely to change the energy 
estimate because of the weak dependence on r\. 

Besides the strength of the magnetic field, there is an uncertainty about the source 
size and the assumption that I ~ r cut . The extent of the region occupied by the relativistic 
particles can be estimated from diffusion length scales. Observations of radio halos in clusters 
such as Coma suggest that the diffusion coefficient of GeV particles in a magnetic field 
of B ~ 2 fiG is D ~ (1-4) x 10 29 cm 2 s" 1 (Schlickeiser et al 1987). The corresponding 
diffusion length scale in a Gyr, assuming a diffusion coefficient of ~ 2 x 10 29 cm 2 s" 1 is 
r D ~ a/ (QDt) ~ 60 (t/10 9 yr) 1 / 2 kpc. We therefore find that t\d ~ r cut ~ 60 kpc for the 
groups in our sample if the groups are 1 Gyr old. For a time scale of t ~ 4Gyr, the particles 
would diffuse out to a region of extent ~ 2r cut . 

As mentioned earlier, our observations are sensitive only up to an angular size of ~ 
40', which corresponds to a length scale (at the redshift of our sample clusters) of 4r cnt . 
If the extent of the region occupied by the relativistic particles is greater than 4r CMt the 
current observations would be insensitive to detect radio emission from such a region. Our 
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observations can, therefore, put a limit on the value of / > 2r cut , and in turn, put a limit on 
the age of the relativistic particles as being larger than ~ 4 Gyr. Our observation therefore 
puts a limit on the epoch of any AGN activity as being at least 4 Gyr before the present 
epoch. Incidentally, a recent proposal of heating the ICM with cosmic rays from AGNs at 
z > 2 (which are also associated with the production of Lithium-6 in halo metal poor stars) 
is consistent with this conclusion (Nath, Madau, Silk 2006). 

This time scale can also be compared with the epoch of formation of structure with 
masses comparable to the total masses of the groups considered here (~ 10 13 M ) in the 
hierarchical structure formation scenario. One can estimate the probability distribution of 
the redshift of formation, defining it to be the epoch by which a fraction of 3/4 of the total 
virial mass has been assembled from merger, given the redshift of observation, that is, the 
present epoch. In the ACDM cosmology, the duration available for evolution after most of 
the mass has been assembled for a 10 13 M is estimated to be however ~4x 10 9 yr (Nath 
2004). We can also find from Figure 1 of Voit et al (2003) that for such a structure, 3/4th of 
its total mass was assembled at ~ t/to ~ 0.7, where to ~ 13.47 Gyr is the present age of the 
universe (Spergel et al 2003); the time available for further evolution is therefore ~ 4 Gyr. 

Our observations therefore put stringent limits on the model of AGN heating in clusters 
for enhancing the entropy of the ICM. The null observations in the case of four groups with 
excess entropy indicate that if AGN heating is the cause of heating the IGM, then it occurred 
well before the mass assembly of these clusters — in other words, the ICM was pre-heated 
before the formation of the clusters. 



4. CONCLUSIONS 

The galaxy groups NGC 524, 720, 3607, and 4697 do not show any evidence of either 
current, or past AGN activities to account for their decreased X-ray luminosities correspond- 
ing to their temperatures. Any possible AGN contribution to the heating of the ICM is a 
factor of ~ 100 lower than that required to explain the decreased X-ray luminosities of these 
groups. Alternatively, any possible AGN activity must have occurred more than ~ 4 Gyr 
ago for the current observations to have not detected its radio emission. The latter scenario 
implies that the ICM was pre-heated before the formation of the clusters. 

We thank Dr. Mitchell C. Begelman for his comments on the manuscript. 
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Table 1. GMRT Observations 



Source 


R.A. (J2000) 


Dec. (J2000) 


RMS 


^cut 




hh mm ss 


o ' " 


(mJy beam -1 ) 


(kpc) 


NGC 524 


01 24 47.7 


+09 32 20 


3.8 


56 


NGC 720 


01 53 00.0 


-13 44 19 


2.1 


65 


NGC 3607 


11 16 54.6 


+ 18 03 07 


1.9 


62 


NGC 4697 


12 48 35.9 


-05 48 03 


1.8 


53 



Note. — Each of the sources were observed at 235 and 610 MHz 
for a total of 6 h. For all the observations the synthesized beam is 
20" X 20" at 235 MHz and 16" X 16" at 610 MHz. The RMS quoted 
is at 235 MHz. The RMS at 610 MHz is ~ 0.2 mJy/beam in all 
cases. 
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log(T) (keV) 



Fig. 1. — X-ray luminosities of poor clusters and groups are shown against their respective 
gas temperatures in keV. The data are from Osmond & Ponman (2004). The V symbols 
are radio-quiet groups and filled circles are radio-loud groups as in Croston et al 2004. The 
straight line is an extrapolation of the L x - T relation from rich clusters, as in, for example, 
Novicki et al 2002. Note that even the radio-quiet groups are below that estimated from this 
extrapolation. 
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Fig. 2. — NGC 4697 : a small section of the radio continuum image at 235 MHz. The 
contour levels are at -4, -2, 2, 4, 6, 8, 10, 15, 20, 30, 40 and 50 times 2 mJy/beam. The 
synthesized beam is 18" X 18". The RMS on this image is 1.8 mJy/beam. The peak flux 
density in this field of view is 1.3 Jy/beam. 
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RIGHT ASCENSION (J2000) 



Fig. 3. — NGC 4697 : same region as in Fig. 2, but at 610 MHz. The contour levels are at 
-4, -2, 2, 4, 6, 8, 10, 15, 20, 30, 40 and 50 times 0.2 mJy/beam. The synthesized beam is 
17" X 13" (PA = -11°). The RMS on this image is 0.15 mJy/beam. The peak flux density 
in this field of view is 0.23 Jy/beam. 



